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ABSTRACT 



A rational method is presented to predict the penetration 
of free-falling objects into deep-sea sediments by combining 
proven empirical theories from the field of soil mechanics with 
known hydrodynamic phenomena. The impact velocity of the ob- 
ject and the shear strength profile and density of the sedi- 
ment are assumed to be known. The penetration problem was 
solved through the use of a computer by equating the work 
done during penetration to the energy of the object falling 
through air and impacting onto a modeled deep-sea sediment. 

The objects were simple geometric shapes ranging in weight 
from 500 to more than 1,000 pounds. The impact velocities 
ranged from zero to twenty feet-per-second . The results are 
compared with full scale tests and recommendations are made 
to extend the method to a water-sediment interface. The 
method successfully predicts the penetration of objects into 
weak, saturated, sediments within the accuracy of the state- 
of-the-art techniques for measuring the sediment mechanical 
properties. The impact duration time was observed to be 
relatively constant and independent of object velocity, shape, 
and weight implying that it may be a unique property of the 
dynamic behavior of a sediment type. 
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INTRODUCTION 



I. 

A. PURTOSE OF INVESTIGATION 

Various types of objects will be lost which fall through 
the water column and impact on the ocean bottom as man con- 
tinues to move more freely on the surface and through the 
oceans. Such objects can range from small and very expen- 
sive instrument packages to fleet submarines. It is important 
to know exactly how far they have penetrated into the bottom 
during search operations since complete penetration below the 
surface of the sediments may require a search by other than 
visual or acoustic methods. 

The velocity and free fall attitude of an object that has 
been located and photographed on the bottom may be estimated 
from the penetration distance if the physical and mechanical 
properties of the sediment are known. The impact velocity 
and attitude could be potentially useful in reconstructing 
the event and also for determining the amount of energy trans- 
ferred to the sediment and radiated from the point of impact. 

The recovery method and force necessary to break a lost 
object free of the bottom during a salvage operation will 
depend to some degree upon how far the object has penetrated 
into the bottom sediments. 

The fields of soil mechanics and foundation engineering 
have been concerned mainly with static situations including 
the prediction of the very slow settlement that a structure 
undergoes as it is erected. The weight is thus applied to 
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the soil in progressively greater amounts. Deep ocean 
bottom installations may be preassembled, floated to their 
site, and allowed to free-fall to a position on the bottom. 
The free-fall method of implantment is potentially useful 
for ocean structures since it decouples the structure from 
the influence of surface motion and tethering lines. Pre- 
diction of an accurate initial penetration depth upon con- 
tact, therefore, becomes a critical design factor. 

An indication of the sediment's mechanical and physical 
properties can be obtained by observing probes of known 
physical configuration as they penetrate the bottom since 
penetration depth depends on sediment strength. 

B. SCOPE OF INVESTIGATION 

The purpose of this investigation was to develop a method 
to predict the penetration depth of objects dropped into a 
simulated deep-sea sediment. Experimental data was available 
to test and verify the analytical method that was developed. 
The full scale experiment employed large, heavy objects and 
it was necessary to let the objects fall through air to 
allow them to reach a velocity at impact equal to the highest 
expected terminal velocity of free-fall through sea water for 
very dense objects. In addition, dropping the objects in 
air eliminated free fall stability problems as the impact 
attitude was effectively controlled and the effects of 
angular momentum on the impact were minimized. The sediment 
used in the experiment was of re-worked continental origin 
characterized by low shear strength and high water content 

typical of deep-sea and pelagic deposits. 
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The initial boundaries and limitations that were imposed 
on the study were designed to simplify the complexities of 
the problem to permit the mechanics to be understood, and 
to retain the ultimate objective of practical application of 
the method at sea. The conditions are summarized below. 

1. The impact velocity at the top of the sediment column 
was in the range from zero to 20 feet-per-second . Zero 
velocity at the top of the sediment column involved releasing 
the object at the sediment surface and its acceleration 
through the soft upper layers before coming to rest. The 
upper limit of 20 feet per second was believed to be a real- 
istic terminal velocity in sea water for most objects. 

2. Free fall was assumed, or gravity was the only down- 
ward force acting on the object. 

3. The object shapes were limited to simple geometries 
with established hydrodynamic properties from theoretical 
calculations or empirical tests. Configurations designed to 
penetrate, such as a coring tool with a long barrel and large 
mass on one end were not considered. 

4. The sediments were taken to be cohesive, with no 
angle of internal friction. They were also assumed to be 
low strength, less than four pounds per square inch, and to 
have high water contents typical of marine sediments. 

5. The study did not include the effect's of any settle- 
ment after initial penetration. 

6. The minimum factors necessary to predict penetration 

depth were considered to be: the impact velocity; geometry 

and mass of the object; and the shear strength profile and 
density of the sediment. 
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C. HISTORICAL BACKGROUND 

The loss of the U.S.S. THRESHER and the problems involved 
in the location of the hull which had penetrated deeply into 
very soft sediments served to emphasize the fact that little 
was known about physical phenomena associated with deep-sea 
sediments and the mechanism of penetration. This disaster 
prompted the original full scale experiments carried out by 
the Naval Civil Engineering Laboratory and used to verify 
the method developed in this study. The THRESHER disaster 
also pointed out the difficulties involved in sampling sedi- 
ments at deep ocean water depths and using only a few samples 
tested in a laboratory to predict as precisely as possible 
the sediment's dynamic behavior in-situ. This study is only 
a step toward understanding all the variables and their 
relationship to each other. 

The following is a review of the state-of-knowledge before 
the investigation and the technical background. The variables 
involved and the development of other penetration equations 
are discussed. 

1 . State-of-Know ledge 

The problem of earth penetration by projectiles is a 
classical one of terminal ballistics and an immense amount of 
both theoretical and experimental literature exists on the 
subject. Empirical attempts to determine the thickness of 
earth embankments needed to protect soldiers were made as 
early as 1742 [Young, 1969]. In recent years ' the free-fall 
penetration problem has been studied to provide protection 
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against artillery, to estimate forces acting on a structure 
during earth impact such as would be required by orbiting 
nuclear power plants that would have to remain intact during 
re-entry and landing on earth, and to determine physical soil 
properties remotely by an acceleration- time record transmitted 
from a projectile penetrometer. The studies have involved 
the fields of engineering, mechanics, fluid mechanics, geology, 
soil mechanics, aerodynamics, wave propagation, as well as 
other disciplines. They were concerned with the problem of 
a body falling with high velocity through air and impacting 
on hard continental soils. Bibliographies of prior investi- 
gations are available [Schmid, 1969 and Young, 1969]. Dynamic 
studies in the field of soil mechanics on land are still in 
their infancy. Whether they can be extended to structures 
and the penetration problems in the deep oceans probably will 
remain a question for some time into the future for four basic 
reasons . 

First, the velocity of an object free-falling through 
water will be much less than the velocity of a similar object 
in air. At very low velocities, from zero to 10 feet per 
second, the inertia of the object may not be important com- 
pared to other factors. In water the added mass of entrained 
fluid must be accounted for, but probably would be ignored in 
air . 

Second, continental soils are quite different from 
deep-sea sediments. Sediments in the ocean basins typically 
have very small grain sizes, often clay sizes of .005 mm and 
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smaller. The rate of sedimentation is measured fractions of 



millimeters per year with very old layered material found 
relatively close to the top of the sediment column. Deep- 
sea sediments are completely saturated with interstititial 
pore water often accounting for half their bulk volume. The 
pore water is not the pure water commonly found on land, but 
the chemically complicated sea water which causes flocculation 
of some sediment particles and the presence of authigenic 
minerals near the sea-sediment interface. Lower in the sedi- 
ment column, diagensis and cementation between particles 
occur forming a loose fabric often with the porosity being 
relatively constant with depth. The effects of the low 2 
degree centigrade temperature and high 0,000 pounds-per- 
square-inch pressure encountered at the median ocean depth 
create an entirely different environment than found on the 
continents. The degree of reworking by benthic organisms 
may also be important in determining the mechanical character- 
istics of the upper layers of the sediments. 

Third, when measuring the mechanical properties of 
the sediments, the amount of disturbance that the sample 
undergoes as a result of the sampling process , handling and 
preparation for testing, temperature and pressure changes, 
and organic growth or decomposition is difficult to precisely 
estimate and may vary from one sampling and testing technique 
to another [Richards, 1961]. 

Fourth, if the sediments in the upper layers are 
treated as a suspension of particles in a fluid, which may 
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be a valid approach to sediment dynamics if the water con- 
tent is over 100 percent, then another problem arises. The 
fluid dynamics of non-newtonian materials have been researched 
for conduit flows such as slurries in pipelines, but very few 
investigations of external flow about a body [Pazwash and 
Robertson, 1969], and unsteady flow at a boundary have been 
made . 

2 . Prior Investigations 

There are only a few records of objects that have 
been lost and recovered. The instances where things have 
been recovered or photographed to document and determine the 
penetration depth together with adequate data about the 
mechanical properties of the sediment are virtually non- 
existent. The same lack of data exists in reports written 
about the few cases where penetration was investigated. The 
reports typically set forth equations with empirical constants 
and contain no information about the controlling sediment 
parameters . 

Much literature and empirical information is available 
about pile driving. The forcing of a long piling into soft 
sediments is, however , quite different from the free-fall 
penetration problem. 

D. TECHNICAL BACKGROUND 

All of the variables were studied, then the development 
of the previous penetration equations was researched before 
an attempt was made to solve the penetration problem. A 
review of the known variables for which data exists was made 
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to develop a method that would be of use at the present 
time using available data and techniques. 

1 . The Variab3.es Involved 

Penetration may be a function of more than 20 differ- 
ent variables and may be described by the following equation: 

f (V , M, A, z , e , r, S, E , v , p T,b, w, s, £ , n 

S (1) 
^ c i c , h / 7 / b i p t a o , h) o 

O WO 

where , 



V Q = impact velocity 
M = mass of the object 

A " cross sectional area of the object 

z = depth of penetration 

r = surface roughness of the object 

S " object shape 

e = void ratio of the sediment 

E = Young's modulus of the sediment 

^ = Poisson's ratio of the sediment 

p = density of the sediment 
s - 1 

T ~ time for penetration to occur 

b = mineralogic and chemical composition of 
the sediment 

p g = viscosity of the sediment 
c = shear strength of the sediment 
k ~ permeability of the sediment 
s ~ sensitivity of the sediment 
w = water content of the sediment 

7 = dilation of the sediment 

8 = grain size of the sediment 
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£ - post depositional changes of the sediment 

p ~ density of the water 
m = viscosity of the water 
h = water depth 

a g = adhesion of the sediment to the object 

Such a list is self explanatory and is not exclusive. 
Some of the above properties are related to others, and the 
parameters may thus be reduced in number by a dimensional 
analysis [Liu, 1969], Some of the properties, however, are 
not amenable to a mathematical analysis, for example the 
grain size distribution and the grain's morphology and 
orientation in the sediment fabric. A phenomological 
approach to the problem may, therefore, be required [Muga , 
1966] . To put all of the variables into a system of 
equations may be an impossible task, and for this reason 
the prediction of penetration into deep-sea sediments may 
remain somewhat of an art and not susceptible to an exact 
solution for some years. 

The deep-sea sediment particles accumulated over 
millions of years occur in many types which reflect unique 
depositional conditions of the earth's history. The 
mechanical properties vary widely, as subsequent to 
deposition the sediments may have undergone secondary 
cementation, diagenesis, recrystallization, and a degree of 
reworking by animals in the top important layers [Smith, 

1969] . Some problems involved in learning more about deep- 
sea sediments are set forth in the following paragraphs. 
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The details of the sea-sediment interface and upper 
sediment layers in the depths of the ocean are only grossly 
known the world over. In the most studied ocean, the North 
Atlantic, Keller [1967], estimates one core has been analyzed 
for sediment mechanical properties for every 30,000 square 
miles of ocean floor. Less data is available in the Pacific. 
The variation of the ocean bottom’s nature and the high cost 
of about 1,000 dollars per foot of core recovery and analysis 
will limit any great expansion of knowledge. 

Much input to the design of foundations on land is 
based on past local experience incorporated into local 
building codes. On land, tests can be conducted and equations 
developed that fit special circumstances in certain areas at 
considerably less expense than in the oceans. Also a past 
history exists for land areas of interest, whereas compara- 
tively little has been known of deep-sea sediments until very 
recent times. Although the horizontal area variability of 
deep-sea sediments may be less than that of continental 
sediments, to empirically categorize the seventy percent of 
the earth's surface covered by water and anticipate a return 
in less than the number of years that have elapsed since 
structures have been built on land would be folly. 

2 . The Development of Previous Penetration Equations 

A review of the development of penetration equations 
emphasizes the importance of the various constants that act 
to control the penetration depth. The discussion that fol- 
lows combines the outlines given by Schmid [1969] , and 
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Young [1969]. First and most basic from Newton's second 
law , 

F ( t ) = JMx(t) (2) 

it can be shown that penetration depth is determined if 
velocity is known as a function of time. In general this is 
not the case for unsteady motion at a boundary. 

Other classical equations are also based on Newton's 
second law where resistance to penetration is a function of 
the instantaneous velocity, V; 

-Mx = C 1 + f 1 (V) + f 2 (V 2 ) + ... f n (V n ) (3) 

which postulates that during the first phase of penetration, 
at high velocity, the motion is governed by resistive forces 
proportional to velocity to some power, and at the final 
stages of penetration the resistance is independent of velo- 
city. If the functions, f ^ , f^, ••• ,f , are taken to be 
constants and the highest power of n assumed to be 2, then 
equation (3) simplifies to; 

-Mx “ q- + / 3 V + 7 V 2 (4) 

where at high velocities the forces controlling penetration 

2 

are proportional to V or analogous to drag force in fluid 
flow. At moderate velocities resistance will be proportional 
to the first power of V, analogous to viscous resistance, and 
at the last stage, proportional to a constant force indepen- 
dent of velocity. According to how the constants in equation 
(4) are treated, the following classical equations can be 
derived after integrating twice with the initial conditions, 
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V ” V Q at t - o, and final conditions, V ~ o at t - T. 



2 

M V 

Robins-Euler z = -2- (constant resisting force) (5) 

2 a 



Poncelet 



Resal 



z 



z 



M 

2 7 



In (1 + 



7 V 



o, 



a 



(resisting forces (6) 
constant and similar 
to fluid drag) 



.. 7 V 

M , , , , 1 o 

— In (1 + 

7 a 



(resisting forces 
similar to fluid 
viscosity and 
fluid drag) 



( 7 ) 



The question relative to the use of the above equations rests 
with the constants, a, ft, y, which are determined by the geo- 
metry and mass of the object and the type of medium being 
penetrated. Schmid [1969] discussed several empirical 
penetration equations from the literature which are vari- 
ations of the above classical equations. The constants were 
determined from high velocity impacts on continental soils. 

Thomason and others [1968] used Poncelet ' s equation 
to determine the amount of charge required for an explosive 
anchor system. 



z 



Pj_ 

2 gib 



In 





( 8 ) 



where 



p = W - weight: of projectile 
A normal frontal area 

a = constant related to the shatter strength 
of the medium (for "mud" a ~ 320 psi) 

b = an inertial coefficient for the material 

i = a form factor, usually 1.0 
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Evaluation of a and b is a problem and unf ortunately 
a comparison of predicted and actual penetration was not 
given. However, the penetration depth's dependence on high 
velocity was illustrated. 

Young [1969] presented a semi-empirical solution to 
the penetration problem derived from 160 full scale tests on 
continental soils with impact velocities ranging from 100 to 
1,000 f eet-per-second . Penetration was treated as three 
phases: (1) impact, (2) moving through the soil, and (3) 

coming to rest. The phenomena of entry is not the same as 
moving through the soil, but the surface effect of entry was 
found to be negligible after the object had penetrated three 
body diameters. The basic form of Young's equation was: 

z = f 1 (S) f 2 (A) f 3 (W) f 4 (V Q ) f 5 (N) (9) 

where 

S ~ a shape parameter of the object 
A ~ cross sectional area of the object 
W “ weight of the object 
V Q ~ impact velocity 

N = constant dependent on soil properties 
averaged over penetration distance 

Two equations were determined for two ranges of velocity, 

one for less than 200 feet per second and one for greater 

than 200 feet per second. The low velocity equation was: 

z = 0.53 (N) (S) (|)ln (1 + 2 x 10" 5 V^) (10) 

and nomograms were given for easy solution. All Young's data 
fitted within + 20% deviation. Again it should be noted that 
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the soil property constant, N, could vary by a factor from 1 
to 50 and "was arbitrarily assigned by those knowledgeable 
in the field". 

If penetration is treated as being resisted by viscous 
forces only, the following results according to Schmid [1969]; 

Mx = A V (t ) (11) 

which has as its solution 



where 



x 



V 

o 

a 



(i 




Am 

M~ 



(12) 



(13) 



When t goes to infinity, the maximum penetration becomes: 

V 

z = — °- (14) 

a 

but again resolving the viscous parameter of the sediment, ju , 
and establishing a realistic time, t, become difficult. 

Soils do not behave as isotropic, linearly elastic 
materials, except possibly at very small strains [D'Appolonia 
and Lambe , 1970]. Penetration into a material causes the 
elastic limit of the material to be exceeded and theories of 
elasticity cannot be applied directly to the problem. 

Theories of elasticity and plasticity applied to problems 
of deformation assume idealized stress strain laws which are 
not well suited to soils and more basically ignore the 
important dimension of time, and assume that the state of 
strain is instantaneous. 

Various viscoelastic models, Figure 1, with combina- 
tions of Kelvin and Maxwell models, of a sediment can be 



21 



\ 




assembled and solved by Laplace or Fourier transforms, if 
time dependence is linear. 



no permanent 
deformation 



KELVIN 

Figure I. KELVIN AND MAXWELL VISCOELASTIC MODELS 

However the problem of determining the necessary soil param- 
eters is non-trivial, and to get a unique solution, several 
sets of parameters are required [Schmid, 1966]. 

Modern advances in rheology and the continued develop- 
ment of the finite element method with appropriate yield 
criteria such as developed by Noh and Wilkins in 1964, and 
outlined by Schmid [3969] , and the analysis of the initial 
settlement problem by D'Appolonia and Larnbe , [1970] , may be 

useful but again the problem of relating all parameters and 
determining the stress-strain behavior and yield criteria of 
a rheologically complicated material such as deep-sea sedi- 
ments, may remain elusive for the near future. 

In very weak, high water content, deep-ocean sedi- 
ments, the penetration problem is somewhat similar to the 
problem of ship slamming. Studies of ship slamming and 
unsteady hydrodynamic impacts at a water surface have mainly 
been concerned with peak pressure, velocity- time relation- 
ships as a function of the form of a structure as it penetrates 



elastic rebound 
permanent deformation 
deformation 



MAXWELL 
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a fluid boundary and not penetration depth. Extensive 
literature and model test records exist [Ochi , 1962]. 

A method to estimate the impact velocity from the 
observed penetration depth was developed by Mandl and Givens 
[1964] after an aircraft was lost and crashed into soft clays 
in northern Canada. The method developed and explicity 
solved an equation of motion for a sphere and cylinder 
penetrating at high velocities into a Leda clay 

- (m ' + M) x = Mg - F - D (15) 

where 

F = the static strength of the clay, a function 
of penetration depth only 

D = a drag force considered a function of 
penetration depth, velocity and time 

Experimental tests of high velocity impacts into clay, 

(greater than 800 feet per second) confirmed that the method 
overpredicted penetration depth. The over-prediction was 
expected since viscous dissipation and elastic wave initi- 
ation were ignored . 

After a detailed examination of mathematical methods 
available for analysis of unsteady motions of an object 
crossing an air-water interface, Moran [1965], concluded that 
none of them could yield a uniformly valid approximation to 
the solution of surface crossing and that the best hope for 
obtaining a reliable estimate was with numerical analysis. 

3 . Known Variables for Which Data Exists 

Since all of the methods reviewed rely totally on the 
ability of the selected variables and constants to correctly 
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parameterize the penetration problem, it was necessary to 
review what is known about the variables. The physical 
properties of the object are assumed to be known but the 
description of the sediment strength is subject to the 
interpretation of empirical tests or, in some instances, 
phenomenological observations. Even though deep-sea sedi- 
ments are very different from continental soils, a review of 
what has been learned about the mechanical behavior of fine- 
grained clays on land and what is known about deep-sea sedi- 
ments is necessary to assemble a workable solution. Also 
the problems of sampling and testing at sea cannot be 
ignored . 

Information about the following parameters; density, 
water content, void ratio, shear strength, and sensitivity, 
which describe the mechanical behavior of deep-sea sediments, 
either exist in the literature or may be easily determined 
through laboratory tests which have been standarized and 
soon may be routinely performed at sea on oceanographic 
vessels . 

The bulk density of deep-sea sediments is determined 
from routine measurement of the weight of a known volume of 
sediment. Typical values range from 74 to 125 pounds per 
cubic foot with the most common sediments in the range from 
78 to 109 pounds per cubic foot [Keller, 1967]. The bulk 
densities of North Pacific sediments are generally lower 
than those of the North Atlantic. 

The void ratio and porosity are a measure of the 
volume of sediment in a sample compared to the volume of 
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voids, which are filled with sea water in the case of deep- 
sea sediments. They are also a measure of particle point to 
point contacts and bonds [Schmid, 1969]. The measurement is 
easily made with a pycnometer that measures that volume of a 
known weight of oven-dried sediment. 

Sensitivity, the ratio of the original undisturbed 
shear strength to a remolded and disturbed shear strength, 
is important to the penetration problem since the object 
moves through disturbed sediment during penetration. Its 
importance is hard to estimate since most samples are dis- 
turbed to some degree. Mitchell and Houston [1969] indicate" 
that in the range of sensitivity from 1 to 8 , typical of 
deep-sea sediments, the mechanism that causes sensitivity is 
chemical cementation, which is not well understood by soil 
mechanists but obviously important. 

For a cohesive sediment, the strength is a measure of 
the sediment's ability to resist mechanical deformation. 

Large masses of cohesive soils generally fail along a sur- 
face and hence the strength of a cohesive soil is defined as 
shear strength. For purely cohesive sediments undergoing 
large deformations, particle to particle interactions which 
would cause an "angle of internal friction" and increase in 
the strength of the sediment are generally assumed not to 
exist. However, as more sophisticated testing techniques 
are developed, this may prove to be an erroneous assumption. 

Typical values of shear strengths were given by Keller 
[1967] to range from .25 to 2.5 pounds-per-square-inch , with 
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the most common range between .5 and 1.5 pounds-per-square- 
inch . The North Atlantic basin sediments have slightly 
higher strengths than the North Pacific. The presence of 
carbonates will increase shear strength, but few carbonates 
are found below 12,000 feet. In the Gulf of Mexico, Bryant 
and others [1967], found shear strengths of silty clay 
turbidites of the abyssal plains to range from .3 to 1.1 
pounds-per-square-inch with the pelagic sediments of the 
continental slopes having strengths up to 5.9 pounds-per- 
square-inch . 

Richards [1967] predicted that a reduction of in-situ 
strength would be caused by sampling and laboratory testing 
and that the reduction might be proportional to the grain 
size, the amount of gases in solution, and temperature in- 
creases that cause increased bacterial action. However, 
samples tested in-situ at 300 meters with a vane device 
correlated well with tests from cores taken from the same 
area a year earlier and tested with a fall cone and laboratory 
vane [Richards, 1967]. Recent experience indicates this may 
have been an isolated instance of agreement. 

Laboratory tests to determine the effects of tempera- 
ture and pressure on shear strength and consolidation have 
been carried out with greatly disturbed, remolded samples 
[Vey and Nelson, 1966]. The results were somewhat incon- 
clusive due to the small number of samples and difficulties 
in instrumentation. The results indicated that there was a 
decrease in direct shear strength with increased environmental 
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pressure for a fine grained sediment with high void ratio. 
The vane shear tests with increased environmental pressure 
showed an ■'^'■fiase in shear strength for plastic sediments 
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be physically contained during _ ting. This pre- 

cludes the use of an unconfined compression test or any test 
that requires trimming the sample before testing. The vane 
shear and fall cone tests can be done on a sample in its 
original container and the vane shear has the added advan- 
tage of easy adaptation to in-situ testing. Kravitz [1970] 
compared the results of Wykham-Farrance vane, Torvane, and 
fall cone tests on weak sediments and found that the results 
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from the fall cone were more repeatable than the Wykham- 
Farrance vane or the Torvane which showed least repeatability. 
These results may have been expected since an important 
factor in vane shear testing is the rate of loading and the 
separation of static and dynamic resistance [Housel, 1959]. 

The vane test's advantages are: First, it can be 

conducted in-situ and a strength-depth profile is easily 
obtained where all other methods require reentry of a hole, 
which may be very difficult in the deep ocean, or discrete 
sampling which disturbs the sediments. Second, it is a 
simple inexpensive technique when compared to sampling and 
laboratory testing. Third, it does not require a system of 
weights or free-falling objects to cause the soil to fail 
and therefore has the potential to be used at sea on a 
rolling, pitching ship. 

The vane test limitations can be summarized as : 

First, it can only be used in uniform, saturated, cohesive 
sediment. Second, when done in-situ it does not obtain a 
sample for evaluation. Third, it imposes a failure surface 
which may not be the weakest surface or relevant to the 
problem being studied. 

Eden [1966] found that the vane will yield higher 
strengths with more consistent results in sensitive clay 
with shear strengths less than ten pounds per square inch 
when compared to a fall cone. This is in conflict with the 
more recent tests of Kravitz and emphasizes the problems in 
testing. Kravitz made statistically significant tests on a 



laboratory prepared montmorillonite sample and Eden tested 
natural clays. 

Shear strength values from a vane were found to be 
15% higher than those obtained by an unconfined compression 
test on Gulf of Mexico slope sediments [Morelock, 1969]. 

This may be partially balanced by the findings of Inderbitzen 
[1970] who found that laboratory vane measurements were 22% 
lower than in-situ measurements made from DEEP QUEST at a 
depth of 3,500 feet. This is in rough agreement with the 
estimation made by Crisp [1968] , of a decrease in shear 
strength due to an expansion of pore water caused by tempera- 
ture and pressure changes. Crisp estimated that an expansion 
of pore water by two percent would occur by raising a sample 
from 12,000 feet deep and 2 degree centigrade temperature to 
surface temperature and pressure. This expansion could 
significantly alter the strength properties of a sediment 
and Crisp showed shear strength should decrease 10 to 15 
percent . 

The in-situ tests conducted by Inderbitzen showed a 
uniform increase in strength with depth in the sediment, but 
the laboratory tests were erratic, suggesting that the differ- 
ence may have been caused by disturbance during sampling and 
preparation for testing. 

Penetrometers offer the great advantage of potentially 
being self-contained instruments. They either record or 
transmit a force versus time curve from which depth of 
penetration by integration and soil properties can be 
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inferred. In order to extrapolate to different shapes and 
velocities, however, the curve's shape must be a character- 
istic of the sediment parameters. There are real problems 
with instrumentation and directional stabilization as the 
penetrometer penetrates since centric impact rarely occurs 
with other than symmetrical objects [Schmid, 1966]. Three 
other problems also complicate the penetrometer approach to 
determining sediment properties in the oceans. First, to 
get deep information, the penetrometer must have enough 
energy to penetrate, and at low terminal velocities in water 
it will probably require a smooth, known, driving force. 
Second, viscous forces may predominate on small penetrometers 
making extrapolation to larger objects very difficult. Third, 
penetrometers do not retrieve samples for analysis of other 
sediment parameters. 

Acoustic sounding and seismic methods do not discrim- 
inate enough of the structure of the upper layers of sedi- 
ments important to the penetration problem to allow their use 
at the present time. In the future the use of high frequency 
pulsed sources mounted on deep vehicles, hold promise for 
better resolution of the upper layers. 

The methods of determining sediment strength that 
have been discussed were developed for use in classical 
foundation engineering where a safety factor based on 
engineering judgement would be applied to the final result. 

The penetration problem is quite different as it requires an 
exact answer since an error factor of two could mean the 
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difference between an object being completely buried below 
the surface or having one half of the object exposed and 
easily located. 



II. METHOD OF ANALYSIS 



Prior to a discussion of the method used to solve the 
penetration problem in the air-to- sediment case and the 
modifications required for use at sea, it is suitable to 
quote directly from Terzaghi and Peck [1948]. 

"Because of the unavoidable undertainties involved 
in the fundamental assumptions of theories and in the 
numerical values of the soil constants, simplicity is 
of much greater importance than accuracy. If a theory 
. is simple, one can readily judge the practical con- 
sequence of various conceivable deviations from the 
assumptions and act accordingly. If a theory is com- 
plicated , it serves no practical use until the results 
are condensed into graphs or tables that permit rapid 
evaluation of the final equations on the basis of 
several different assumptions." 

A. ASSUMPTIONS THAT PERTAIN TO THE SEDIMENT 

Deep-sea sediments are not well defined in the literature 
dealing with soil mechanics. Therefore a list of important 
assumptions considered in developing the penetration pre- 
diction method are stated : 

1. The sediment is completely saturated with sea water. 

2. The water content of the sediment is 75 percent or 
more, as is typical of marine deposits in the deep ocean 
basins. 
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3. Organic material that would tend to bind the sediment 
particles together is either absent or is present in very 
small quantities. 

4. The sediment will fail along its weakest plane. 

5. The only shear resistance is cohesion. 

6. The initial penetration takes place in such a short 
time that the effect of permeability is nil. In a study of 
time-dependent deformation of clays Barden [1969] noted that 
the effect of adsorbed water on individual particles was that 
pore pressure equals the applied pressure at the first 
application of load. 

7. Any lithif ication or cementation that may have developed 
along individual particle-to-part icle contacts is broken at the 
time of impact. This implies a portion of the sediment is 
remolded at impact. If the particles do not react with one 
another by other than mechanical means, the mixture will 
behave as a viscous fluid [Peck and others, 1953, Pazwash and 
Robertson, 1969] . 

8. Shear strength varies in an irregular manner with 
depth. The variation of velocity gradients is greater in the 
deep ocean basins than over shallow areas [Ewing and Nafe, 

1963] which implies that deep ocean sediments have been 
deposited in layers. This is confirmed by visual observation 
of deep-sea core samples which often show horizontal layers 
of different colors and textures. When studying in-situ 
strengths, Inderbitzen [1970] found that horizontal vari- 
ations were much less than vertical variations. He also 



33 



found that the percentage of variation in shear strength 
decreased with depth , which may indicate that biological and 
chemical processes are taking place in the upper layers and 
minimal consolidation occurs. 

9. The bulk wet density of the sediment can be considered 
constant with depth in the surficial sediments. The shear 
strength of a saturated sediment depends to a large extent on 
the structural arrangement of the particles and Scott [1967] 
points out that soils can have different structures at the 
same density, but may not have the same structure at differ- 
ent densities, which implies that shear strength is not a 
unique function of density. X-ray diffraction and polarized 
light microscope techniques cannot reveal the details of 
particle configuration or orientation in fine-grained marine 
sediments. However, electron microphotographs of sediments 
from the Gulf of Mexico show them characterized by a loose, 
open arrangement of randomly oriented particles and that 
particle rearrangement resulting from gravitational consoli- 
dation is not apparent [Bowles, 1969]. 

10. The energy transmitted to the sediments in the form 
of elastic or gravity waves from the initial impact shock is 
negligible in comparison to the total energy dissipated during 
penetration. Impacts on soils and fluid media have been 
studied for a long time, but there is a very difficult instru- 
mentation problem with detecting energy transmitted away from 
the impact site. According to Hunter [in Schmid, 1966] the 
Hertz impact theory, which refers to impacts in the elastic 



34 



range, states that the portion of total energy that is 
vibrational energy is very small as long as the impact 
velocity is less than the velocity of sound in the medium. 
Schmid [1966] measured the distant dynamic: disturbances 
caused by a small penetrometer impacting into materials whose 
rheological properties were known or easily tested, such as, 
greases, wax, and homogeneous soils. High frequency vibration 
was detected in the hard materials and it was determined by 
use of surface accelerometers that the largest amount of 
energy propagated from the point of impact was 4.3 percent 
of the total energy available. His calculated propagation 
speed suggested that the motions measured at the surface 
were Stoneley waves. Ochi and Bledsoe [1962] when conducting 
model tests of ship hull forms slamming onto a water surface 
found that forces were very large for .01 seconds and pressure 
due to impact phenomena were over within .07 seconds, yet the 
model's velocity decrease during this time was less than 4 
percent. Chuang [1970] in a similar investigation, did not 
detect any significant acoustic pressure variations at a 
distance away from the impact point and concluded that water 
could be considered incompressible for impact velocities from 
0 to 50 feet per second. 

B. METHOD OF SOLUTION FOR THE AIR-TO-MUD CASE 

A rational, phenomenological approach was used to solve 
the penetration problem rather than an empirical method 
because of the poor state of knowledge that exists relative 
to actual behavior of deep-sea sediments. V7ith engineering 
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judgment, this method may be extended to different objects 
and different cases more easily than an empirical approach to 
the problem. 

1 . Role of Hydrodynamics 

The sediment was treated as a homogeneous, incompressible 
fluid to enable an estimate to be made of (1) the momentum 
transferred from the body to the sediment, and (2) the drag 
forces exerted on the object at high velocities. If the sedi- 
ment was a true fluid it would have an infinite strain for an 
applied stress. Therefore, the forces required to initiate 
failure in the sediment and bring the body to rest were taken 
as the forces calculated using the ultimate bearing capacity 
formula developed by Prandtl and Terzaghi. 

Treating deep-sea sediments as fluids has been an 
approach applied to the study of turbidity flows and slumping 
on the continental slopes. The term "spontaneous liquefaction" 
is often used when referring to saturated clays and implies 
fluid properties. Briefly, JMorgenstern [1967] states that 
liquefaction occurs when the sediment is loaded to the collapse 
stress of the sediment fabric and the load is then transferred 
to the interstitial pore water. The increase in pore pressure 
then produces hydraulic gradients which further disrupt the 
sediment fabric causing the solids and pore water to flow 
macroscopically as a viscous fluid. Houston and Mitchell 
[1969] describe the fluid phenomena as occurring in sensitive 
clays which contract during shear because of their open 
fabric. If the clays are saturated, the contraction tendency 
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results in a large transfer of normal stress from the clay 
skeleton to the pore water causing pore pressure to increase. 

In saturated deep-sea sediments this transfer of normal stress 
may be almost instantaneous. Continuing the description of 
Houston and Mitchell [1969] , as the remoulding caused by con- 
tinuing strain becomes more pronounced, the effective normal 
stress may fall to a very low value. Very little strain is 
required to break bonds at particle contacts and there is 
some build-up of pore pressure even before a significant 
number of bonds are broken because of slight elastic compression 
of the fabric and incompressibility of the pore water. The 
resistance of sensitive clays due to dilational energy was 
determined to be small or absent by Houston and Mitchell [1969] , 
and Walker [1969]. Because, first, at high void ratios little 
or no tendency toward dilation was observed during shear, even 
at high strains, and second, the effective stress had been so 
reduced at high strains the mobilized resistance of the clay 
fabric was small. Mandl and Givens [1964] when conducting 
high velocity penetration tests into a clay, judged that the 
compressibility phase of impact lasted until the surface 
boundary breaks and that the static yield strength of the 
clay was exceeded by many orders of magnitude under impact 
and the clays behaved in a dynamic manner similar to a fluid. 
Schmid [1966] when conducting low velocity free-fall impact 
experiments on a clay also observed that liquefaction occurred. 

Another reason to treat deep-sea sediments as a fluid 
was that their bulk wet density has a narrow range of variation 
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with depth in the upper layers of sediment and the assumption 
of a median density over the range of expected penetration is 
a good approximation. 

2 . Ultimate Bearing Capacity Equation 

The bearing capacity of a soil is a measure of its 
ability to support a load relative to its limiting equilibrium. 
According to Harr [1966] , the first rational approach to 
developing a method to predict ultimate load carrying capacity 
of a soil was provided by Terzaghi who recognized that if the 
three dominant factors contributing to bearing capacity, the 
weight of the soil, the effects of surcharge, and the strength 
parameter of the soil, could be assessed, a good estimate of 
bearing capacity would be obtained. Terzaghi chose as a model 
of the failure mechanism of the soil the solution to the prob- 
lem of penetration of circular punches into metal as developed 
by Prandtl and Reissner. In the original solution Prandtl 
was able to neglect the effect of the weight of the metal 
moved to the influence of surface loadings on the surrounding 
material. When the material was assumed to be isotropic and 
homogeneous with no angle of internal friction, Prandtl 's 
solution yielded, 

q d = > D t + 5 * 14C (16) 

On the basis of experiments and Prandtl* s solution, Terzaghi 
and Peck [1948] derived the following semi-empirical equation 
for the bearing capacity of a circular footing. 

= 1.3c N N + .6 7 N t 

c f q 7 



q 



(17) 



where q = pressure at surface to cause failure (psi) 

c = cohesion or shear strength of the soil (psi) 

y ~ unit weight of the soil 

= depth of footing below the surface of the soil (in) 

N^,N ,N = factors which depend on the angle of internal 
^ Y friction and depth of footing 

The first terra represents a contribution due to the 
shear strength of the soil, the second the effect of material 
removed by putting the footing below the surface (equivalent 
to buoyancy in water) , and the third term, the effect of an 
angle of internal friction (assumed to be zero for deep-sea 
sediments) . The equation is modified slightly for square or 
oblong footings and Meyerhof [1961] gives constants for wedge 
and cone shaped footings. More complete discussions of the 
bearing capacity equation and slight variations in the values 
of the constants can be found in texts on soil mechanics 
[Terzaghi and Peck, 1943; Harr, 1966; and Hough, 1969]. 
Terzaghi 1 s equation was the closest to Prandtl's (which 
requires theoretical failure surfaces) and it has been 
empirically verified for many years. It was therefore 
selected as a reliable relationship for estimating the static 
strength of deep-sea sediments. The values of the constants 
and the shapes used in equation (17) are given in Table I. 

In the bearing capacity equation, the shear resist- 
ance mobilized on the failure surface of any surcharge 
material is neglected, hence the equation is conservative 
and relates only to events occurring below the footing or 
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penetrating object. buring penetration, static resistance 
will occur in the undisturbed portion of the sediment below 
the deepest penetrating part, the other surfaces of the 
penetrating object are subjected to remolded sediment or 
sediment already set in motion by the penetrating object. 
Therefore the bearing capacity equation constants used in 
penetration prediction were considered to be for shallow 
footings for all depths of penetration. 

3 . The Character of Failure 

Before a final method of solution was selected , 
motion pictures of penetrations into soft sediments and still 
photographs of soil failures were studied. In the motion 
pictures, the sediment was observed to move laterally outward 
from the object and have somewhat of a separation point from 
the object. While dropping plates vertically into a simu- 
lated sediment, Erchul [1S68] found that polymer coatings to 
reduce drag on the plates had little effect at high entrance 
velocities and postulated that the sediment may not have been 
in contact with the plates. He also observed that at low 
velocities the lubricants did reduce adhesion to the plate 
and considerably greater penetration occurred with lubricated 
plates compared to unlubricated ones. 

To find a qualitative indication of the effect of 
fluid in interstitial spaces, Rowe and others [1962] dropped 
probes into finely ground (less than 40 microns) silica dust. 
At atmospheric pressure the dust displaces laterally and 
craters were formed. In a vacuum with no air in the inter- 
stitial spaces, resistance to penetration was 100 percent 
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greater and the dust did not move laterally and no mounds 
or craters were formed. 



The motion pictures of penetrating objects in the 



test pit revealed little splatter and little surface dis- 
turbance at a distance from the penetration point. Clark 
and Robertson [1965] observed motion pictures of a buoyant 
body exiting from water and came to the conclusion that 
energy lost due to surface effects was small. Moran [1965] 
also arrived at the same conclusion. 

4 . Equation of Motion 

Since penetration is a dynamic phenomena involving 
the unsteady motion of a body, the first attempt at evalu- 
ating the penetration depth involved writing an equation of 
motion for the object as it penetrated the sediment. The 
final solution required a work energy approach that allowed 
the change in velocity and penetration depth to be determined 
incrementally with a computerized numerical scheme without 
the requirement to integrate the equation of motion and 
establish the instantaneous relationship between force, mass 
and acceleration. 



characterized as a fluid with a static strength are shown in 
Figure 2. The equation of motion can be written as follows: 



Where the applied forces are: weight of the object, the sedi- 

ment static strength, the buoyancy provided by the sediment, 
the force to overcome the inertia of the sediment, and the 
drag forces. 



The forces acting on an object penetrating a sediment 



Ma 




Applied Forces 



(18) 
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Figure 2. FORCES ACTING ON A PENETRATING OBJECT 

The change of the body's motion as it penetrates the 
sediment is the equivalent of the resistance encountered by 
the body, or the energy involved in the change is equal to 
that expended in propelling the body. Momentum is trans- 
ferred to a surrounding fluid by; (1) a certain volume of 
fluid being accelerated uniformly, (2) the fluid being put 
into turbulent motion, (3) compression waves, (4) gravity 
waves or spray [Hoerner , 1958]. All of the momentum is 
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eventually consumed by viscous friction. The last two modes 
of momentum transfer have been assumed to be negligible for 
the velocity ranges and object sizes under consideration and 
were discussed previously. 

The treatment of the drag force, or the momentum put 
into the fluid by turbulent motion and viscous drag, outlined 
below, follows that given by Mandl and Givens [1964], 

The drag forces cannot be evaluated easily for unsteady 
motion through a fluid boundary. The coefficient of drag can 
be approximated by; 



1 3 

C = (— + C ) 
d l Re c dJ • 



(19) 



The first term represents the contribution due to viscous 

drag; the second is that due to form drag; 0 is a constant 

depending on the shape of the body and is usually represented 

by the first term in an infinite series (Oseen's approximation). 

64 

For a circular disk placed broadside to a stream 0 ~ — , while 

7T 

for a sphere 0 = 24 (Lamb) . The term C. corresponds to the 

■Mto 

asymptotic case as Reynolds number approaches infinity in 
Navier-Stokes flows, it depends on geometry and not the fluid 
dynamic properties. Inspection of experimentally obtained 
drag coefficients versus Reynolds number relationships for 
bluff bodies show that is a slowly varying function at 
large Reynolds numbers. For bluff bodies with sharp edges 
nearly all the drag is form drag. The experimentally deter- 
mined values of C, in the case of bodies with sharp edges is 

^ DO 

approximently 1.0. For bluff bodies without sharp edges, 
such as spheres, the form drag also depends on the position 
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of flow separation which in turn depends on the critical 
range of Reynolds numbers. The equation of motion then 
becomes : 

(m’+M)x = Mg - p Q - 7 g A x - A ~ Ps AC^x 2 (20) 

If the viscosity of clay is taken to be on the order 

of 3 x 10 , lbs - sec/in [Whitman in Mandl and Givens, 1964] , 

the Reynolds number for the objects that were under study with 

a three foot maximum dimension and with a velocity of ten 

2 

feet-per-sccond , is on the order of 2 x 10 . Admittedly a 
constant value of viscosity for a rheologically complicated 
substance such as a deep-sea clay is a gross over simplifi- 
cation, but the order of magnitude of the terms indicates that 
the viscous drag contribution to the equation of motion may 
be ignored [Mandl and Givens, 1964] . Also, towards the 
terminal phase of motion where the effect of viscosity is 
liable to increase rapidly with decrease in the rate of 
strain, the influence of viscosity will be negligible since 
the velocities are small and the other terms in the equation 
are larger by several orders of magnitude.- Since form drag 
represents pressure forces on the body which are significant 
at the higher velocities, it could not be neglected. 

5 . Added Mass Effects 

Another reason to treat the sediment as a fluid was 
to have a method available to estimate the volume and velo- 
city of sediment displaced as the object penetrated. Womack 
[1967] verified that the amount of saturated clay moved 
during penetration was equal to the size of hole created. 
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The energy required to move the sediment is conserved since 
once the sediment is put in motion by the object it will 
continue to move away from the object creating a larger 
crater as the object is slowed by other forces. It is very 
important though , because the drag term of the equation of 
motion is dependent on velocity squared and the inertial 
effect of the sediment motion will alter the velocity- time 
history of penetration which has significant effects depend- 
ing upon the shape and impact velocity of the object. 

Added mass is caused by hydrodynamic pressure forces 
which are proportional to the instantaneous values of accel- 
eration and according to Wendel [1950] appear to have the 
same value in a viscous or inviscid fluid as long as small 
bodies are not involved. The added masses were determined 
according to equations and methods given by Saunders [1957] , 
which are tabulated in Table II. 

The assumptions pertaining to the fluid approximation 
of deep-sea sediments can be summarized as follows and arc 
the same as given by Moran [1965] in a study of water entry 
and exit. 

1. Irrotational flow - viscous effects on the pres- 
sure felt by the body are generally small and the variations 
in inviscid forces are more important than viscous drag. 

2. Incompressible flow - impact velocities are well 
below the speed of sound. 

3. Surface effects are negligible. 



45 



The problem with determining depth of penetration 
from equation (20) is that it cannot be integrated during 
the unsteady motion of an object transiting a fluid boundary 
from one fluid to another. This is true for the following 
reasons: (1) m' , the added mass and mass of the displaced 

sediment, is a function of depth of penetration for an object 
whose submerged portion changes shape as it submerges; (2) 
the dominant term, that of static soil strength, varies with 
the depth in a highly irregular manner; (3) the location of 
the air-to-sediment boundary condition is unknown a priori; 
and finally, (4) the flow is unsteady. 

6 . Transformation of the Equation to Work-Energy 

The equation of motion may easily be l'esolved into 
an equation of energy and work which is time independent and 
a function of increments of distance. The development is 
given below. 




d 

dt 



(21) 



Multiplication by dx gives 



d , dx , 
dt dt ' 



dx 



,dx , 



2> x ax - h ax * m a q£) 



( 22 ) 



The first integral with respect to displacement is : 



J F x dx = 1/2 M (|^) + C (23) 

The left hand term represents the work done on M during the 
interval dx and can be evaluated if is a known function 
of x. The right hand term represents the corresponding 
change in kinetic energy of the object. 
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Kinetic energy is a scalar quantity that depends only 
on the mass and magnitude of the velocity. Since the velocity 
is squared it is always positive. The total mechanical energy 
of the body is equal to: 

E = KE + PE (24) 



For an object penetrating the sediment and coming to rest, 
the energy of the object can be accounted for by: 



E 1 + e in e out 



(25) 



This simple balance states that the total energy E^ of the 
object at the beginning of an interval plus any energy put 
into the object minus any energy taken out must equal the 
final energy, < °£ the object at the end of the interval. 
When the object has lost its energy in overcoming buoyancy, 
static soil strength, inertia of the sediment, and drag, it 
will stop its motion. 

A particular reference level for the energies is not 
needed since energy is relative. The potential energy 
depends on an arbitrary selection of datum for zero potential. 
Also the kinetic energy is arbitrarily expressed relative to 
the condition at zero velocity. Hence any physical measure- 
ment over an interval describes a change in energy. It 
follows that reference to an arbitrary datum used for the 
expression of energy will cancel. 

The work-energy method approach to solving a problem 
is powerful but dangerous because an answer is always obtained 
with no indication that the answer is correct. The method 
balances work to equal energy available and if a portion of 
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of the energy available is not considered , or incorrectly 
computed, work is still done and an answer results. In the 
unsteady case of penetration, using an iterative solution 
almost requires a computer because errors in computation will 
not be indicated but will appear as energy gained or lost. 

The answer will not appear wrong unless a gross error was 
made . 

In analysis of the penetration problem only mechani- 
cal energy was considered, with the sediment and object 
being a conservative system in the gravitational field. 
Thermal, chemical, electrical and atomic energies were 
neglected . 

7 . Method of Solution Used with a Computer 

The numerical iteration scheme used to solve the 
penetration problem was the same for each object shape. 
Briefly, it divided the object into sections or slices in 
the same plane as the sediment surface. Then the work 
necessary to submerge each section and the work done by 
the sum of the submerged sections, the forces due to the 
inertia of added mass and displaced mass of the sediment, 
the static strength of the sediment, and drag in moving an 
incremental distance equal to a section thickness were 
calculated and subtracted from the kinetic energy and change 
in potential energy of the object. Penetration depth was 
determined when the total kinetic energy of the object 
equaled zero. 

A description in mathematical terms follows and a 
block outline of the computer solution is shown in Figure 5. 
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Constants and necessary approximations for each object are 
listed in Table I. The computer programs are included as 
Appendixes C through H. 

Refer to Figure 3 for an explanation of the 
definitions that follow. 




d 






Position 1 




Position 2 



Figure 3. DEFINITIONS 

v. = average velocity of the object as it transits 
distance x from position 0 to position 1. 
(Assumes a linear velocity gradient) 



49 



v 2 “ average velocity of the object as it transits 
distance x from position 1 to position 2 

M = mass of the object 

m^ ~ added mass at section 1, plus the mass of 
displaced fluid at section 1 

mi, ~ added mass at section 2, plus the mass of 
displaced fluid at section 2 

A m = m^ - m^ 

U = work done 

Ufo = work done against buoyancy forces as object 
transits from position 1 to position 2 

U = work done against sediment static strength as 
the objects transits from position 1 to 
position 2 

Ihj = work done against drag forces as the object 
transits from position 1 to position 2 

AU = the change in energy of the body as it transits 
from position 1 to position 2, not including 
change in potential energy 

= total energy of the object at position 1 

= total energy of the object at position 2 

A PE = change in potential energy of the object over 



distance x 



From equation (22) 



v 



2 




m ' v d v 



(26) 



v 



1 



U~ 1/2 m ' v 2 - 1/2 m 1 v ^ 



(27) 



it follows that 



AU = 1/2 m^ v2 - 1/2 m^ v^ + 1/2 m'v^ + u b + u s + U d 



(28) 
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NOTE: 1. Am' was accelerated from a static condition, zero 

velocity to the velocity at the end of the section 
over which it was computed. 

2. Am' was assumed to be zero after the object had 
penetrated to its maximum dimension since 
separation of the sediment from the object was 
apparent at this point. 

The incremental distance, x, is fixed, and the various masses 
are known as is the work necessary to move the object from 
position one to position two. Velocity at position two is 
unknown. Therefore from the work energy relationship 

E 2 = E 1 + PE - AU = 1/2 Mv 2 (29) 

Combining equation (28) with (23) and solving for gives 



2 (APE + E 1 + 1/2 m 



v. 



1 



V- 



- U, - U - 



b 



V 



(M + Am' + 






(30) 



In the iteration scheme, when" v^ equals zero, the sum of the 
increments, x, is the penetration depth. 

The work done against buoyancy forces is a straight- 
forward calculation that is discussed in elementary fluid 
mechanics texts. In the computer iteration scheme where the 
sections of an object are taken to be slices with vertical 
sides it becomes; 

2 

~ 1/2 7 A x (for the section being submerged) 

2 

plus 7 A x (for all sections below the one 

being submerged) 

The work done against drag forces was taken to be: 

= 1/2 pCj Av^x (31) 

where 7 “ bulk wet density of the sediment 

C d = drag coefficient for the approximate object 
shape that is below the sediment surface 



A = the maximum cross sectional area of the object 
below the sediment surface and normal to the 
velocity vector 

v ~ average velocity over the distance moved, 
assuming a linear variation in velocity. 

The work done against the sediment static strength 
was calculated from the ultimate bearing capacity equation. 

U s = c f N c A c x (32) 

where C_^ = a constant depending on footing shape 

N ^ = a constant depending on the type of sediment 
and footing profile 

c = shear strength of the sediment 

When the penetrating object had a flat plane forcing 
into the sediment such as i\ cylinder impacting on end, the 
cross sectional area. A, of the flat area and the shear 
strength, c, of the sediment at the penetration level of the 
flat area were used in the equation above. In the case of 
objects with rounded or sloping sides, a judgment had to be 
made of a cross sectional area and representative shear 
strength. The area selected was the maximum cross sectional 
area normal to the velocity vector at the surface of the 
sediment, or the maximum cross sectional area of the object 
if the area had penetrated beneath the sediment surface. 

This method of approximating a sloping or curved surface by 
a projected horizontal area appears valid for blunt objects 
[Meyerhof, 1967]. Since the shear strength characteristics 
of deep-sea sediments in the upper layers appear to change 
radically with depth in the sediment column, a shear strength 
for equation (32) was determined by selecting an average 
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depth according to the zones of plastic failure implied by 
equations (16) and (17). This was selected as one-half of 
the penetration depth until the maximum cross sectional area 
had been reached, then it became the depth of penetration 
minus one half the penetration depth when the maximum cross 
sectional area was coincident with the surface of the sedi- 
ment as illustrated in figure (4) . 




Figure 4 . DETERMINING THE DEPTH OF "c" 

;>3 



FOR NON-FLAT OBJECTS 



The method described was adequate for curved and 
conical objects and approximates a representative strength 
in the theoretical zones of plastic failure in the sediment. 

An alternative and more refined method of selecting A and c, 
would have been to use the differences in the cross sectional 
area from section to section and apply the value of c at the 
depth of each section. This increased the complexity of the 
iteration procedure , hence was attempted only when there 
were discontinuities in the penetrating object's shape such 
as might occur with a truncated cone. 

The distance of each iteration was chosen only to 
conserve computer time and approximate the expected accuracy 
of the field tests. Larger and smaller iteration intervals 
were tested for effects on the fina] depth and as long as 
the interval was less than one-half of the distance between 
shear strength measurements, convergence was always apparent. 
With a linear shear strength profile, the interval is limited 
only by the desired accuracy of the solution. 

The equations for added mass are given in Table II. 

A flow chart of the basic program is given in Figure 5. 

Appendix A presents a set of curves which qualitatively 
illustrate the relative influence of each force in the pene- 
tration equation. They were obtained by using the sediment 
shear strength profile typical of the world's deep-ocean 
basins as determined by Iloag [1970] and shown in Figure 13. 
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Table I. CONSTANTS AND SHAPES USED FOR COMPUTATIONS 



